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ABSTRACT 


This  report  describes  a  technique  for  calculating  the  input  impedance  of  ideal, 
infinitely  long  TEM  horn  antennas  used  to  radiate  ultrawideband  signals.  A  simple 
potential  calculation  is  used  to  determine  the  fields  of  the  TEM  mode  of  the  anterma, 
from  which  the  input  impedance  is  derived.  This  result  is  useful  for  matching  TEM 
horns  to  the  output  impedance  of  the  excitation  source  to  avoid  undesirable  reflected 
signals. 
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Input  Impedance  Calculation  for 
TEM  Horn  Antennas 


EXECUTIVE  SUMMARY 

A  simple  numerical  technique  is  described  for  calculating  the  input  impedance  of 
TEM  horn  antennas  for  ultrawideband  (UWB)  applications.  UWB  antennas  are 
required  to  radiate  subnanosecond  pulses  with  minimal  distortion.  The  required 
bandwidth  is  from  a  few  hundred  megahertz  up  to  10  gigahertz.  TEM  horn  antennas 
have  emerged  as  potential  candidates  for  directional  UWB  applications  such  as 
impulse  radar.  The  input  impedance  of  the  horn  is  typically  matched  to  the  output 
impedance  of  the  source.  At  the  DSTO  Wideband  Test  Facility,  the  antenna  geometry 
has  previously  been  varied  empirically  to  produce  the  required  input  impedance  due 
to  the  unreliability  of  alternative  theoretical  techniques  based  on  cordormal  mapping 
which  was  not  corrected  until  recently.  This  simple  analysis  circumvents  this  time 
consuming  process  and  is  readily  implemented. 
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1.  Introduction 

Recently  there  has  been  a  growing  interest  in  the  possible  defence  applications  of 
ultrawideband  (UWB)  radiation.  These  applications  are  currently  being  investigated 
at  the  DSTO  Wideband  Test  Facility.  The  generation  and  transmission  of  these 
ultrawideband  emissions  such  as  impulse  radiation,  requires  a  fast  electrical  switch 
to  generate  a  subnanosecond  transient  which  directly  excites  a  wideband  antenna.  If 
the  excitation  pulse,  or  more  commonly  its  time  derivative,  is  to  be  faithfully  radiated 
in  the  far  field  zone,  then  a  nondispersive  antenna  is  required.  This  last  requirement  is 
a  significant  restriction,  and  rules  out  log-periodic  antennas  commonly  used  as 
compact,  wideband  antennas. 

One  class  of  antenna  that  is  both  wideband  and  nondispersive  is  the  TEM  horn,  as 
shown  in  Figure  1.1.  Two  quasi-triangular  plates  are  connected  together  at  one  vertex 
to  form  the  feed.  A  spherical  TEM  wave  propagates  along  the  antenna  structure, 
confined  mainly  between  the  plates.  When  the  TEM  mode  reaches  the  end  of  the 
antenna  structure,  conversion  to  higher  order  propagating  spherical  TM  and  TE 
modes  occurs  because  the  TEM  mode  can  no  longer  be  supported.  This  type  of 
antenna  has  been  shown  in  the  past  to  faithfully  transmit  the  time  derivative  of  the 
incident  pulse  in  the  far  field  zone  [1]. 


Figure  1.1:  Conventional  TEM  horn  geometry 


The  input  impedance  of  TEM  horns  is  dependent  both  on  the  flare  angle  0a  and  the 
apex  angle  of  the  plates,  (|)a.  For  this  analysis,  the  length  of  the  horn  is  assumed  to  be 
much  longer  than  the  spatial  width  of  the  excitation  pulse,  so  that  reflections  from  the 
end  of  the  horn  can  be  ignored.  For  the  100  ps  pulses  normally  used  in  the  DSTO 
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Wideband  Test  Facility,  the  spatial  pulse  width  is  3  cm.  This  is  small  compared  to 
the  TEM  horn  antenna  plates  which  are  usually  1  m  long.  Therefore,  the  input 
impedance  of  the  horn  is  the  same  as  for  an  ideal  infinitely  long  antenna  for 
frequencies  of  practical  interest.  Note  that  the  impedance  of  an  antenna  which 
supports  only  the  TEM  mode  is  independent  of  frequency.  The  input  impedance  of 
the  anterma  is  generally  required  to  match  that  of  the  voltage  source  (typically  50  Q). 
Previously,  this  match  has  been  obtained  empirically.  This  report  describes  a  simple 
analysis  based  on  a  finite  difference  approach  to  calculate  the  TEM  mode  electric  field 
of  an  ideal  infinitely  long  horn.  From  the  calculated  field  distribution  for  given  angles 
0A  and  (|)A/  the  impedance  is  calculated  from  the  power  transmitted  through  the 
cmterma. 

Usually,  <1)A  is  given  and  0a  is  to  be  found  based  on  the  desired  anterma  impedance. 
Note  that  as  (j)A  increases  the  beamwidth  in  the  horizontal  plane  increases.  To 
maintain  the  matching  to  the  source  impedance,  0a  must  increase.  This  has  the  effect 
of  decreasing  the  beamwidth  in  the  vertical  plane  as  <|)a  increases.  Antenna 
beamwidths  are  not  predicted  by  this  analysis,  but  the  ability  to  calculate  the  value  of 
0A  required  to  match  the  antenna  to  the  source  impedance  for  a  given  <j)A  removes 
one  degree  of  freedom  from  the  horn  design.  This  considerably  simplifies  the  design 
process. 

Carrel  [2]  was  the  first  to  develop  an  analysis  intended  to  calculate  the  input 
impedance  of  infinitely  long  TEM  horn  antennas.  This  technique  relied  upon 
conformal  transformations  to  yield  a  simple  analytic  form  for  the  impedance. 
Unfortunately,  Carrel's  analysis  has  been  found  to  be  incorrect  due  to  an  error  in 
applying  the  conformal  transformations.  Lambert  et  al  [3]  have  recently  produced  a 
revised  version  of  Carrel's  analysis  which  hc»s  been  shown  to  yield  the  correct  results. 
Note  that  until  Lambert  et  al  developed  a  corrected  conformal  mapping  analysis,  this 
technique  was  considered  too  unreliable  to  use.  This  report  presents  a  simple  two- 
dimensional  analysis  based  on  a  finite-difference  technique  which  does  not  require 
conformal  mapping  and  is  particularly  simple  to  implement.  In  addition,  antennas 
with  arbitrary  cross-sections  in  the  (0,  (j))  plane  can  be  analysed  since  the  technique  is 
two-dimensional.  This  work  arose  from  the  application  of  a  finite-difference  time- 
domain  analysis  for  ultra  wideband  antenna  modelling  [4].  The  results  of  this  finite- 
difference  technique  are  in  good  agreement  with  experimentally  measured  input 
impedance  values  for  several  TEM  horns. 
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2.  Theory 


2.1  The  Finite  Difference  Calculation 

From  Collin  [5],  a  TEM  mode  electric  field  in  spherical  geometry  can  be  related  to 
scalar  functions  g  and  O  by 


where  O(0,<t))  satisfies 


=  0 


and  g(r)  follows 


^  +  =  0 


with  k  =  (oVlxoEo-  Therefore 


g(r)  =  e-f^ 


0^^  +  sin 


in0cos0^  +  ^]<J>(0/<|))  =  0 


Note  that  Equation  1  does  not  contain  any  dependence  on  k  and  hence  frequency. 
Therefore,  this  analysis  applies  directly  in  either  the  frequency  or  time  domain.  The 
potential  O(0,(|))  can  be  represented  as  0(i,y)  on  a  spherical  finite  difference  mesh 
discretised  in  0  and  (j)  using 

7C  ,  K  . 

0i  =  2  ■  *  ^0  =  2  '^‘1* 

where 


Ae  = 


2(N0+1) 


and  A,|,  = 


for  i  =  0,  1,.,.,  Ne,  and  j  =  0,  1,...,  N(|,.  Symmetry  allows  the  entire  antenna  to  be 
analysed  by  considering  only  the  quadrant  0  =  0  to  7c/2  and  ^  =  n/2to  37r/2.  The 
spherical  coordinate  system  is  shown  m  Figure  2.1. 
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Figure  2.1:  Coordinate  system  for  TEM  horn  input  impedance  analysis 


The  antenna  plate  was  assumed  to  have  zero  thickness  and  extend  over  (j)  from  j  =  0 
to  j  =  jpiate  the  theta  value  corresponding  to  i  =  ipiate-  This  requirement  for 
discretisation  for  the  finite  difference  representation  of  Equation  1  meant  that  the  flare 
and  apex  plate  angles  of  the  antenna,  0a  and  (t)A/  were  rounded  off  to  the  nearest 
(6,-,(t)y)  values  on  the  mesh.  The  effect  of  this  rounding  could  be  minimised  by  a 
judicious  choice  of  Ng  and  N<j,. 

Expressing  Equation  1  in  its  finite  difference  form  leads  to 


sin  0 

sin  0  cos  0’ 

sin  0 

'sin  0  cos  01  , .  , 

■sin2  0  1  ■ 

Ae 

_  Ae  2  . 

+ 

> 

<X> 

Li,  -  2 

2  2 
^<t). 

+  0{i,j+l)\+^{i,j-l)\ 

A.  A. 


=  0 


This  can  be  written  as 


where 


A<D  =  b 

=  ^{hj) 


(2) 


(3) 
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Ait(k-l) 


Ajt(fc+1) 


1_ 

2 


_  sin  9  r sin  6  cos  9~ 

^k{k+jmax)  ~  Ae  Aq  2 

_  sin  6  r sin  6  cos  6~ 

^k(k-jmax)  ~  Ae  [  A9  '  2 

otherwise 

Aki  =  0 

for  k  =  i  (N(|)+l)  +  j.  For  free  space,  =  0  for  all  k,  however  Ay  and  bj^  will  also 
include  the  boundary  conditions  which  are  described  in  the  next  section. 

Equation  3  can  be  solved  using  successive  overrelaxation  based  on  a  Gauss-Seidel 
iterative  algorithm.  An  initial  estimate  of  the  potential  0(i,;)  is  made  for  all  i,  j.  The 
value  for  the  potential  at  each  grid  point  is  updated  using 


(m+1) 

4 


,  V  (m+1) 
(1-CO)  —  -  CO 


l<k 


(m) 

^kin 


l>k 


(m)  , 

Ay  Ojt  -  bfc 


Afefc 


The  iteration  was  halted  when  the  average  magnitude  of  the  correction  to  the 
potential  over  the  entire  mesh  was  less  than  lO'^  typically.  This  condition  was  used 
to  specify  adequate  convergence. 

The  relaxation  parameter  co  is  given  by 


CO 


2 

1  + 


where  |i  is  the  spectral  radius  of  the  block  Jacobi  matrix  B  =  -D‘l  (L  +  U) 
corresponding  to  A  =  D  +  L  +  U  [6].  D,  L,  and  U  are  the  diagonal,  upper  triangular, 
and  lower  triangular  elements  of  A,  respectively.  The  spectral  radius  p  is  the  largest 
eigenvalue  of  B,  which  can  be  readily  calculated  [7]. 
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Once  calculated,  the  potential  can  be  used  to  determine  the  electric  and  therefore  the 
magnetic  field  of  the  TEM  mode  in  the  antenna. 

2.2  Boundary  Conditions 

Boundary  conditions  were  required  at  the  edges  of  the  mesh.  For  0  -  0  (/  =  0),  the 
potential  O  was  zero,  since  the  anterma  was  represented  by  a  positively  charged  plate 
at  0  =  +0A  and  a  negatively  charged  plate  at  0  =  -  0a-  Therefore  when  i  =  0, 

Afcfc  =1  and  Ak(k-l)  =  Ak(k-l)  =  A.k(k-jmax)  =  ^k{k+jmax)  =  0 

At  0  =  0A  {i  =  ipiate)>  the  potential  was  set  to  Vo/2,  where  Vq  is  the  total  potential 
difference  across  the  plates.  Therefore  for  i  =  ipiater 

Akk  =1  and  Ak(k-l)  =  Ak(k-l)  =  ^k{k-jmax)  =  ^k{k+jmax)  =  0 

and 

bfc  =  Vo/2 

Since  the  TEM  horn  anterma  is  symmetric  about  the  y-z  plane  in  figure  2.1, 

3<I>  K  3k 

^  =  0  at  <|)  =  2  O' =  0)  and  (t)  =  y  O' = 

Hence  which  means  that  for  j  =  0 

1 

Ak(fc-i-i)  =  2 


and  for  j  =  N(|, 

1 

Ak(k-i)  = 


The  final  boundary  condition  was  derived  for  the  potential  0  =  0  position  at  the  top  of 
the  spherical  mesh.  At  this  point,  the  potential  O  =  ^top-  As  0  0,  Equation  1  infers 

that 


32<D 

9(1)2 


0 


If  the  potential  in  the  vicinity  of  0  =  0  (i  =  Ne)  is  represented  in  terms  of  (()  harmonics 
as 
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o 


w  =  0 


then 


-n^  <I)„  0 

This  implies  that  =  0  near  0  =  0  for  n  7^0.  The  n  =  0  component  can  be  found  by 
averaging  the  nearest  neighbour  values.  From  continuity  of  the  electric  field  in  free 
space,  the  n  =  0  component  of  O  at  0  =  0  can  be  approximated  by  the  value  at  0  =  A©, 
since 


dOo 

^0  — ^  0  as  0  — ^  0. 


Therefore,  for  0  =  Aq  (i  =  Ne)  Equation  2  can  be  written  as 


r sin  0  cos  0'! 

■sin2  0  1  ■ 

Lao  -  2  J-2«>(Ne,j) 

2  2 
.  ^<t>- 

+  O(N0,;+1)  \  +  0(N9,;-1)  ^ 

A.  A,, 


_  sin  0  r sin  0  cos  0' 

"  ■  Aq  [  Ae  2  . 


where 


^top  ~  NT.+i  ^(Ne,;) 

;  =  0 


r4 

t 


Therefore,  no  modifications  to  Ajd  need  be  made,  and  b^  is  now  given  by 


bfc  =  - 


sin  0 


Ae 


sin  0  cos  0 
.  Ae  ^  2  . 


N^+1 


I 

i  =  o 


<E>(N9,;) 


for  i  =  Ne. 


2.3  Electric  Field  Calculation 

The  electric  field  quantities  <XQ{i,j)  and  a(^{i,j)  which  are  related  to  the  electric  field 
through 
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Ee= 


and  E((,  =  Vq 


«0(w) 

r 


are  calculated  in  the  centre  of  the  finite  difference  mesh  cells  using  the  values  of  the 
potential  at  the  vertices  of  the  cells  as  follows 

^  +  Mil)  -  ^m,i)  -  0(i+l,/+l) 

MW  -  -  30  2  Ae 

/•  v  1  3<1>  0(i,/)  +  g>(i+l,y)  -  g>(t'3+l)  -  <l>(i+l,/+l) 

®<t>(*d)  030  2  sin  0  A,), 

ae(0,(|))  and  a(|,(0,(j)))  can  be  found  for  the  other  three  quadrants  using 


a<t)(ex/<|);)  =  -  a(^(-0xv4>;)  =  - 

and 

a0(0i,<t)/)  =  ae(-0i,<t>;)  =  =  o.Q{-Qi,-^j) 

The  electric  field  can  now  be  calculated  for  the  TEM  mode  of  the  anteima  from  the 
potential  calculated  on  the  finite  difference  mesh.  The  input  impedance  can  now  be 
calculated  from  the  TEM  mode  field. 


2.4  Input  Impedance  Calculation 

The  input  impedance  was  determined  by  calculating  the  power,  P,  transmitted 
through  a  spherical  surface  enclosing  the  anterma  as  follows 

7C  271 

P  =  J  J  ^  ^ sin  0  dcj)  d0 

0*^  0*^ 

where  |  E  |  is  the  magnitude  of  the  electric  field,  is  the  impedance  of  free  space, 
and  a  is  the  radius  of  the  spherical  surface.  Using  symmetry  and  substituting  for  the 
electric  field,  P  can  be  expressed  as 

2  Tt/2  2>%I2 

Vq  f  r  22 

^  ^  ^  J  J  («0  +  “4,)  sin  0  d<t)  d0 
0  %/2 

The  power  can  also  be  related  to  the  voltage  across  the  plates,  Vq,  and  the  antenna 
input  impedance,  Za,  using 
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Therefore 


Co 


71/2  371/2 


J  1 

0  7t/2 


2  2 

(ttg  +  a^ij)  sin  0  d(|)  d0 


In  most  practical  situations,  the  TEM  horn  will  be  fed  by  50  Q  coax  with  an  applied 
voltage  V,„.  If  the  antenna  is  matched  to  the  feed  cable,  then  Vq  =  V,„.  However, 
when  there  is  a  mismatch,  the  voltage  on  the  antenna  plates  is  related  to  the  applied 
voltage  by 


Vo 


2Za 


Za  +  Z, 


•V,-, 


coax 


where  Zcoax  -  50 

The  finite  difference  potential  calculation  therefore  allows  the  input  impedance  of  the 
TEM  horn  to  be  calculated  from  the  TEM  mode  fields  and  from  power  considerations. 
The  only  approximations  are  due  to  the  finite  difference  mesh  discretisation  and  the 
small  rounding  of  the  antenna  flare  and  apex  angles,  0a  and  <t)A/  as  described  in 
Section  2.1. 


3.  Results 

The  antenna  input  impedance  was  measured  for  plate  apex  angles  of  20°,  30°,  and 
45°  using  a  single  triangular  fin  supported  above  a  ground  plane,  with  the  anterma 
fed  from  underneath  the  ground  plane  via  an  N-type  connector.  The  impedance  of 
this  configuration  is  half  the  impedance  of  a  full  horn.  The  impedance  was 
determined  for  a  range  of  flare  angles,  0a  from  the  reflection  coefficient.  The  flare 
angle  was  limited  to  30°  due  to  mechanical  constraints  on  the  feed.  Reflections  from 
the  aperture  of  the  antenna  were  time-gated  out  so  that  the  measurements  could  be 
related  to  infinitely  long  antennas.  A  Hewlett-Packard  HP54120A  18  GHz  digitising 
oscilloscope  was  used  with  an  HP  54121A  four-channel  test  set  which  possessed  a 
time-domain  reflectometry  (TDR)  capability.  The  reflection  coefficient,  p,  for  each 
anterma  configuration  was  determined  by  measuring  the  amplitude  of  the  reflected 
voltage,  Ef,  and  dividing  it  by  the  height  of  the  incident  voltage  step  E,-,  which  was 
measured  as  200  mV. 
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Er 


The  error  in  p  was  calculated  from  the  error  in  E^.  and  E,-,  (t(E),  using 

1  + 

^(p)  = 

The  error  a(E)  was  set  to  the  cursor  width  of  ±1.56  mV. 

The  antenna  input  impedance  was  calculated  from  the  reflection  coefficient  using 

Za  = 

where  Zq  =  50  Q  is  the  impedance  of  the  oscilloscope  input.  The  impedance  error 
C!{ZX)  was  calculated  from  the  error  in  the  reflection  coefficient 

2Za 

<J(Za)  =  (1  .  p)2  c^(P) 

6Xp  6Xp 

Table  3.1  shows  the  measured  input  impedance  Z^  and  the  error,  c{Z^  )  together 

with  the  theoretical  results  Z^  from  the  finite  difference  potential  calculation,  for 
various  flare  angles,  0A/  for  the  20°  apex  angle  TEM  horn.  Tables  3.2  and  3.3  show 
the  same  parameters  for  30°  and  45°  apex  angle  horns.  These  results  are  plotted  in 
Figures  3.1,  3.2,  and  3.3.  The  theoretical  results  are  the  continuous  curves  and  the 
points  correspond  to  measured  impedance  values.  The  agreement  between  the 
theoretical  and  experimental  results  is  usually  within  experimental  error.  The 
calculated  impedance  was  found  to  be  systematically  greater  than  the  measured 
impedance  for  the  45°  apex  angle  horn,  as  shown  in  Table  3.3  and  Figure  3.3.  This 
was  probably  due  to  the  mechanical  difficulties  involved  in  rigidly  supporting  the 
large  triangular  plate  above  the  ground  plane. 

The  number  of  finite  difference  mesh  points,  Ne  and  N,j„  were  typically  in  the  range 
100-150.  The  tolerance  for  the  overrelaxation  scheme  was  10"^.  Under  these 
conditions,  the  run  time  on  a  Sun  SPARC  5  workstation  was  approximately 
2-5  minutes  for  each  flare  and  apex  angle. 
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Figure  3.3;  Input  impedance  against  flare  angle  6 a  for  TEM  horn  with  45°  apex  angle 


4.  Conclusions 

A  relatively  simple  analysis  has  been  presented  which  is  useful  for  calculating  the 
input  impedance  of  TEM  horn  antennas  suitable  for  radiating  ultrawideband  pulses. 
For  ideal,  infinitely-long  antermas  the  input  impedance  is  independent  of  frequency. 
This  approximation  is  valid  provided  the  length  of  the  horn  is  much  greater  than  the 
spatial  pulse  width.  This  is  certainly  the  case  for  the  TEM  horns  used  in  the  DSTO 
Wideband  Test  Facility  which  are  usually  1  m  in  length.  A  finite  difference 
calculation  in  spherical  geometry  was  used  to  calculate  the  TEM  mode  potential 
which  could  then  be  used  to  calculate  the  input  impedance  of  the  antenna  with 
considerable  accuracy.  This  approach  is  simpler  than  alternative  conformal  mapping 
techniques  [2],  [3]  which  until  recently  have  been  considered  as  unreliable. 
Experimental  data  agrees  well  with  the  theoretical  predictions  for  the  TEM  horn  input 
impedance.  For  plate  apex  angles  of  20°  and  30°,  the  agreement  between  the 
measured  and  predicted  impedance  values  is  generally  better  than  a  few  percent.  The 
discrepancy  between  theoretical  and  experimental  results  was  somwhat  worse  for  the 
horn  with  a  45°  apex  angle  at  approximately  10%.  This  was  attributed  to  mechanical 
difficulties  in  supporting  the  plate.  This  analysis  can  now  be  used  to  simplify  the 
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design  of  TEM  horn  antennas  by  allowing  0a  to  be  calculated  for  a  given  <t)A  and 
excitation  source  output  impedance. 
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